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ABSTRACT: We perform a detailed semi-analytical analysis of the electroweak phase tran- 
sition (EWPT) property in NMSSM, which serves as a good benchmark model in which 
the 126 GeV Higgs mixes with a singlet. In this case, a strongly first order electroweak 
phase transition (SFOEWPT) is achieved by the tree-level effects and the phase transition 
strength ye is determined by the vacua energy gap at T = 0. We make an anatomy of the 
energy gap at both tree-level and loop-level and extract out a dimensionless phase transi- 
tion parameter R,, = 4kv,/A,, which can replace A, in the parameterization and affect the 
light CP odd and even Higgs spectra. We find that SFOEWPT only occurs in Ry, ~ —1 
and positive R, < O(10), which in the non-PQ limit case would prefer either a relatively 
light CP odd or CP even Higgs boson ~ (60, 100) GeV, therefore serves as a smoking gun 
signal and requires new search strategies at the LHC. 
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1 Introduction and motivations 


In the last two years, the ATLAS and CMS collaborations have established the discovery of 
the long expected standard model (SM)-like Higgs boson h, with a significance up to 6.1 and 
6.9 a, respectively [1]. This new resonance has a relatively light mass mp ~ 126 GeV, and its 
observed production or decay rate is close to the SM one. With more data accumulating, 
we would enter into the territory of precise understanding of EWSB mechanism. In an 
orthogonal direction, one may wonder its impacts on the weak scale cosmology, in particular 
the corresponding thermal property: the nature of electroweak phase transition (EWPT). 
It is not only a big question of early cosmology per se, it will also help us understand the 
origin of baryon asymmetry in a sense that a strongly first order EWPT (SFOEWPT) is 
required for successful electroweak baryogenesis (EWBG). 

Baryogenesis has a close relation with the Higgs physics and moreover the successful 
baryogenesis implies a non-standard Higgs boson (For discussions on CP violation, see Ref. 
[2]). Broadly speaking, with the current LHC data on Higgs production and decay, we can 


specify three classes of SFOEWPT models based on their discovery potential through Higgs 
physics. The first class is that there is a colored or electric charged particle which couples 
to the 126 GeV Higgs boson. In this case, new particles which alter the Higgs production 
or decay through gluon fusion or di-photon decay channel will change the Higgs effective 
potential and potentially enhance the EWPT strength [3-11]. Comprehensive studies have 
been carried out after the LHC data and it is found out that for a single particle, a 
SFOEWPT requires enhanced gluon fusion production rate and suppressed Higgs di-photon 
decay width [3, 4]. This problem can be cured if one introduces another particle with its 
loop contributions opposite to the first one while the EWPT strength is enhanced [4]. The 
second class is that we have a singlet scalar which couples to the Higgs but never develops 
a VEV [12]. In this case, future precision electroweak and Higgs measurements would 
constrain the overall kinematical renormalization of the 126 GeV Higgs induced by this 
model. The last class is that the extra scalar gets a VEV and mixes with the Higgs (or 
through a tadpole term which is essentially the same [13]) or there are multi-Higgs [14]. In 
this case, it is the mixing effect that changes the Higgs physics properties. Investigating 
its genuine features clearly is an important task. 

Supersymmetry (SUSY) is a well motivated example among beyond SM models and 
it can also provide the SFOEWPT for successful EWBG. For instance, in the minimal 
supersymmetric SM (MSSM) the significant Higgs-stop coupling can lead to an acceptable 
EWPT strength in a tiny window, given a well organized stop sector [15]. In light of the 
recent LHC Higgs discovery and stop exclusion, this window has been severely constrained 
and essentially ruled out [16-18] (For a remedy, see [4]). In this class of model, one challenge 
after the Higgs discovery is to lift m, with the least fine-tuning while still accommodate 
the Higgs constrains. One simple extension is the NMSSM, which provides a large tree- 
level Higgs mass and a natural solution to the u problem. With an extra singlet in the 
Higgs sector, it is conceivable SFOEWPT is still viable in NMSSM and we are curious 
on the phase transition patterns constrained by the current data. In this article, we have 
studied this problem in great detail and found out a critical parameter R, = 4kv,/A, 
where SFOEWPT only occurs in R, ~ —1 and positive R, < O(10), which in turn would 
prefer a lighter CP odd or even Higgs boson. 

This paper is organized as follows. In Section 2 we review the NMSSM in detail, 
including both its Higgs potential in zero and finite temperature. In Section 3, we first 
analyze the SFOEWPT in NMSSM semi-analytically through both the tree level and loop 
effects, and then provide the numerical results of the parameter scan which includes all the 
current experimental and strong electroweak phase transition conditions. We also show 
the corresponding particle spectra patterns, LHC observations and dark matter in section 
4. Finally we conclude and give a discussion in Section 5 and some necessary details in the 
paper are given in the Appendices. 


2 The NMSSM at Zero and Non-zero Temperature 


As mentioned in the introduction, the NMSSM can accommodate natural SUSY with the 
current data constraints and provides viable dark matter candidates, thus receives a lot of 


attention. It also provides a good benchmark model that the 125 GeV Higgs mixes with 
a singlet getting a VEV , therefore provides tree-level cubic terms in the Higgs effective 
potential to enhance the strength of EWPT. In this section we review the basic formulas 
for the Higgs effective potential setup without and with the finite temperature corrections. 


2.1 Tree-level Higgs Potential 


All of the above eminent features of the NMSSM are traced back to the Higgs sector, which 
in the Z3— invariant form, is written as 


Wz D ASHu - Ha + m 
1 
—Lsoft = (aam G H4S + gh Ans + he.) ‘ (2.1) 


After S getting a VEV vs = (S) around the electroweak (EW) scale, an effective y-term 
is generated. This is the original motivation of the singlet extension. But as a great 
bonus, the model provides a Higgs quartic term at tree level, i.e., \?|H?H9|?, which can 
be significant for a large A ~ 1 and moreover a small tan ~ 1. As a consequence, the 
tree-level mass of the SM-like Higgs boson becomes: 


m = (mz, cos? 23 + d? sin? 28) + ÔMmix, (2.2) 


where 6mmix stands for the mixing effects on Higgs boson mass. It can be sizable (typically 
a few GeV), positive or negative depending on the mass order of the SM-like Higgs boson 
mass among the neutral Higgs bosons [19, 22]. If hsm is the lightest one (Hı —scenario), the 
effect is a reduction. In contrast, if hsm is the next lightest one (H2—scenario), the effect 
is an enhancement. After the LHC Higgs 126 GeV Higgs discovery, those spectra pattern 
of those two scenarios have been studied intensively [20-37]. Here we will investigate the 
status of EWPT in these two scenarios separately. 

For later convenience, we give the complete tree-level Higgs potential, which consists 
of the D-, F- and the soft SUSY breaking terms: 

Vo = |AHy: Ha — «S?|’ + |AS]? (Hi Ha + HLH) 

g + t 2 g t 2 
+ q Hu = H} Ha) =F D HH 
+ mi, H} Ha + mi, H} Hu + m32|s/? 


1 
+ (AA, Hy: HaS + gh AnS +h.c.), (2.3) 


where HI = (H9, H7), HT = (H}, He) and g? = (93 + g?) /2. Here we will not discuss 
the CP violation aspects of the electroweak baryogenesis and assume A, A), K, A, E€ R for 
simplicity. An SU(2)xU(1) gauge is chosen such that at the physical vacuum 


vt = (Ht) =0, vu, = (He) ERT. (2.4) 


v- = (M7) = 0 is a local minimum provided that the charged Higgs bosons have the 
positive mass square. Moreover, \A, > 0 and Ax > 0 are assumed to realize vg, vg € RT. 


The angle 8 is defined as in the MSSM: 
vq = v cos b, vy = vsin p, (2.5) 


with v ~ 174 GeV. 
For the Higgs mass square structure, we decompose the Higgs fields as follows [19]. 


1 
H? =v, + a cos 3 + Sy sin 8) + — (P; cos 8 + G° sin £), 


B Z 


1 
H? =vg + —(— S; sin B + S2 cos B) + Ji 


V2 
S3 + iP2 


V2 ’ 


where G° is the Goldstone boson. In this basis, the doublet block has already been ap- 
proximately diagonalized and S% is the SM-like which carries electroweak VEV among the 
doublets. In the basis (S1, S2,.53), the elements of the CP even Higgs mass squared matrix 
elements (M2);; are given by 


(P, sin 8 — G° cos £), 


S =vs+ (2.6) 


u = M3 + (m2, — dv?) sin? 26, 
1 
12 = = 57 — X°?) sin 48, 


(M5) 
(M5) 
(M2)13 = —(M3 sin 28 + 2AKv?) cos 285; 
(Ms) 
(M5) 


M2)20 = M? cos? 28 + A?v? sin? 26, 
1 
M2 23 = zC Må sin? 28 + Any? — 2Akv? sin 28), 
Us 
2 LP sin? v\? 2,2 le. oa: 
(M3)33 = JMA sinf 28 = + 4k U$ + KÅkUs — JARY sin 28. (2.7) 
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where M% = 2Av,(A) + Kvs)/sin 26 defines the largest scale among these elements and 
is the heavy CP odd Higgs mass. We can introduce an auxiliary parameter C4 = 1 — 
A) sin 23/2 — «sin 26/X to measure the mixing between singlet and doublet, i.e., M3; = 
2C4Apv (The mixing (M2)12 can be safely neglected for moderate large tan 8). The other 
light Higgs diagonal mass square (M2)33 can be written as 


1 s 1 1 
(M)33 = = M4 sin? 28 =) a Any? | 1+ — | — =Anv? sin 28, 
4 Us Rk 2 


1 1 
-3 (M5)23 (=) + 4K? v? (1 ) H dy? — Akv? sin 28 (2.8) 
Us 


and we will use this formula again and again in later discussions. Here Rẹ = 4kv;/A, is a 
dimensionless critical variable defined for SFOEWPT. 


For the CP odd Higgs boson, A, is theoretically upper bounded in order to keep the 
CP odd singlet-like Higgs mass squared (M74,)29 positive [38]: 


(Mb)u1 = MĀ, 
12K?v2 


1 a: 
(Mp)22 = JMA sin? 28 (=) = au sin 26 — Ro 


S 


il 
(M>)12 = 5 M4 sin 26 (=) (2.9) 
2 Us 
where (M?) corresponds to the mass squared M2 of the only CP odd Higgs in the 
MSSM. 


2.2 Effective potential at Finite Temperature 


The starting point for the perturbative analysis of EWPT is the finite temperature effective 
potential. Up to one-loop order, it takes the form of 


V (p, T) = Volpi) + Vily T) + Vaaisy (Yu, T). (2.10) 


where yı, l= d, u, s are the classical field variables corresponding to H9, H?, and S. The 
tree-level part Vo follows directly from the Higgs potential in Eq. (2.3). We realize that 
our analysis at this precision may be subject to corrections from high order and the issue 
of gauge dependence [39, 40] and a more complete analysis is left to a future study. 

The one loop part V; consists of the Coleman-Weinberg potential at zero temperature 
and thermal corrections at finite temperature [41]: 
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—)25i 7. 2 A 
H= > Se mile) (in wee >) + = > nidi (=) (2.11) 


where 7 runs over all particles in the model, with each having degrees of freedom n;, field- 


dependent mass m;(yı) and spin s;. J; is the thermal integral function Jg(p) for bosons 
(fermions) 


Jg rly’) = f dx x? In(1 F eV ®t), (2.12) 
0 


It tends to be zero in the non-relativistic limit, i.e., y2 >> 1. By contrast, it has the high 
temperature expansion and in particular gives rise to the well known thermal cubic term 
in Eq. (2.11), given that 7 is a boson. Here we work in the Landau gauge and in the DR 
scheme. As for Vaaisy, it is the daisy resummation contributions from the longitudinal 
components of gauge bosons and the scalar bosons [42—44] 


T 
Vaaisy = ED" +o n(m (Y1, T) = m3(¢¥1)), (2.13) 
b 


where mọ is the thermal mass. 
Finally, it should be emphasized that in our analysis, we use the three VEVs w as the 
inputs and eliminate the Higgs soft masses via the minimization conditions for the three 


field variables yı. Concretely, at one-loop order, they are given by 


=2 — 
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Quo Oey = 
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(2.14) 


3 Electroweak Phase Transition in the NMSSM 


With previous preparations, in this section we study EWPT in this model. It is well 
known that successful EWBG requires a SFOEWPT, namely ye = ve(Te)/Te Z 0.9. For a 
dedicated study of this condition based on gauge invariant quantities, see Ref. [39]. Here 
Te is the critical temperature of SFOEWPT, with order parameter ve(Te). In the SM, the 
lattice simulation indicates that its EWPT is actually a crossover, failing to achieve any 
jumps in terms of order parameter. In the MSSM, in particular after the discovery of the 
126 GeV Higgs boson, a single light stop alone would be ruled out by the current Higgs 
data because of too large enhancement on the Higgs production rate from gluon fusion. 
Nevertheless, a second colored light scalar can not only reduce the gluon Higgs effective 
operator, but also enhance the EWPT strength [4]. While the NMSSM, by virtue of its 
tree-level effects, provides a simple way to enhance ye. Such effects have been studied by 
many groups before [45-52], but a detailed general analysis of SFOEWPT after the Higgs 
discovery is still absent and we fill the gap in this paper. 

We will first introduce three types of EWPT and then propose a new way to investigate 
ye from the zero temperature Higgs effective potential. Following this way, we make an 
anatomy of each type, giving semi-analytical treatment on tree-level effects and qualitative 
analysis of loop corrections. It is found that the latter plays a robust role in SFOEWPT, 
despite of the dominated tree-level effects. 


3.1 Vacua energy gap and SFOEWPT 


The NMSSM contains three Higgs fields and thus possesses a rich vacua structure, which 
leads to a variety of EWPT patterns. There are mainly three patterns [49, 53], classified 
by the course of the phase transition from the symmetric phase Qo to the EW symmetry 
breaking phase Qgw (here we denote various phases with their VEV’s): 


Type-I: Qo > Qs > Ogw At high temperature, the universe is in the symmetric phase. 
As the universe cools down, it may transit to the vacuum locating at the singlet sub- 
space, i.e., Qs. As T further decreases to the critical temperature Te, Qs degenerates 
with Qpw and then the universe transits into the phase Qpw £. 


‘Notice that if < S >= 0 in Qgw this would induce a symmetry non-restoration effect in the S direction 
[54]. 


Type-IT: Qo => Qy => Qew Type-Il transition passes the intermediate phase Qy with Hu 
developing a VEV first. Here only the first step are relevant for EWPT, which recovers 
the SM case except that S contributes to the thermal cubic terms [41]. However, in 
this case, SFOEWPT requires large interactions between Higgs and singlet which 
induces sizable mixing between the two. This will change the transition type into 
Type-III. Generally, Type-II is hardly strong so we will not discuss this type. 


Type-III: Qo = gw The EW symmetry breaking vacuum develops first, and thus the 
universe in the symmetry phase transits directly into the phase rw. It is worth 
pointing out that even though the transition does not undergo other phases, there 
are still extra local minima at T = 0. In particular, there usually exits local minima 
in the singlet subspace which makes the vacua structure Type-I-like. We will turn 
back to this point in later discussions. 


Vacua structure at T = 0 should encode information on EWPT. For instance, the 
effective potential in Type-I is likely to have a metastable vacuum Qs besides the EW 
vacuum Qpw, with vacua energy gap AV = Vg — Vew. The 7T-dependent terms in the 
finite temperature potential need to smooth out this gap as T increases, until the critical 
temperature. Accordingly, a smaller AV(T = 0) may imply a lower Te thus a larger ye. 
This conjecture is confirmed by our final numerical results shown in Fig. 1. In the three- 
dimension field space, developing an analytical expression for ye is mission impossible, 
except for some simplified cases like in the PQ-limit [50, 51]. Therefore, our observation 
provides an important guideline for achieving a larger yc. 
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Figure 1. Strong correlation between the strength of EWPT and the vacua energy gap AV at 
T = 0. Left panel: Hı— scenario; Right panel: Hə— scenario. The vertical line stands for the 
AV = v?m? /4 = 1.18 x 108 GeV limit. 


This general correlation between the phase transition strength and the vacua energy 
gap AV can be understood in the AV — 0 limit. Let’s consider the vacua energy gap 


between the symmetry phase and the broken phase: 


AV = (Veym — Vew (vo)) |r=0 
> Vsym(Te) = Vew (TL = 0, vo) 
= Vew(Te, ve) — Vew (T = 0, vo) 
OV 

~ Teg T = Try) (3.1) 
where in the second line, we have used the approximation Vsym(Te) = Vsym(T = 0) in the 
AV — 0 limit which holds exactly if the symmetric phase is the origin of @ or is a good 
approximation if the symmetric vacuum is of a weakly coupled singlet. The third line 
comes from the degenerate vacua condition at the critical temperature. In the last line, we 
set Uc ~ vo and use the Lagrange mean value theorem and T, € [0, Te]. Thus, 


Ve OV 1 
a (oF non) AV ee 


The dependence of the phase transition strength on AV can be revealed more explicitly 
in the simplified model of NMSSM in the PQ-limit [50]. In this model, 


vt [z 207m? 
ee eee (33) 


For a extremely small AV, ve © v, and the phase transition strength can be rewritten as: 


o 242m2 4 
Zena aL N (3.4) 
Te (m2 + A202 . %4 )2 2AV 


where E is the coefficient of thermal cubic term ¢°T. 


3.2 Anatomy of Type-I 


In this subsection we will lead the way to generate a smaller energy gap AV (T = 0) in 
Type-I. A full understanding requires analysis at both tree-level and loop-level. We also 
give numerical results, which are consistent with those semi-analytical understandings. 


3.2.1 Tree-level analysis 


First of all, let’s investigate the vacuum energy Vew of Qgw. Substituting Eq. (2.14) into 
Eq. (2.3), one can eliminate the Higgs soft masses in Vew, and then Vgw can be divided 
into three parts, Viy, Vey and VES . The first part denotes the contribution completely 
from the Higgs doublets 

2 20 2 


vom 
Vw = a (g’v” cos? 28 + A7v" sin? 28) ~ — i h, 


(3.5) 


To derive the second approximation we have used nothing but Eq. (2.2) where we neglect 
the mixing effects for the Higgs boson mass. As one can see, Ve is definitely negative. 


Moreover, its value is determined by the Higgs quartic coupling thus related to the SM- 
like Higgs boson mass, mp œ~ 126 GeV. Therefore, this part is almost fixed to be around 
—1.18 x 108 GeV*. The second part Ve is the contribution from the singlet sector, taking 
the form of 1 
Vw = — 3 Ans — Kv. (3.6) 
The third part ye is a result of the doublet-singlet mixing, and it can be casted into a 
simple form, 
VES = —Cap’e. (3.7) 


As we have mentioned before, C4 measures the mixing between singlet and doublet M3; = 
2C'4yv. From the current Higgs data, we expect this auxiliary parameter C’, is usually 
much smaller than 1 [28] 7, so the singlet or the SM-like Higgs mass is not largely pushed 
down in the Hə— or the H,—scenario respectively. This fact will help us to simplify 
discussions and furthermore find out a crucial variable R, which has a close relation with 
the vacua energy gap AV and the EWPT strength ye. 

Next we discuss Vg, the tree-level potential energy of the absolute minimum us in the 
singlet subspace °. To compare with Vw, it is convenient to eliminate m2, in favor of vg 
through the third equation of Eq. (2.14), rewritten as 


me, = —C4 Nv? — KAgUs — 2K7 02. (3.8) 
Then from the potential with only singlet S 
V(S) =m? + Ék AS? + nt, (3.9) 
we can get 
Vg = [—K Ax (vs — 2us/3) — kK? (2v? — u2)| u2 — CaA? uo’. (3.10) 


It is also illustrative to express Vg in terms of the inputs only, 


At 


Vas 
S~ 384K2 


(1+ V1—8a,)° (1+ VI— 82, — li (3.11) 


which holds for 2, = m2/A?2 < 1/8, see more details in Appendix. A. Vs is definitely 
negative for £w < 1/9. Moreover, it is an even function of both A, and «. For a given 
£u, Eq. (3.11) indicates that Vg becomes more negative as A2/« (or A,/K to some extent) 
increases. 

With all the above expressions of potential energy, we proceed to discuss the vacua 
energy gap at tree level, which is found to be related to the deviation of us from vs. To 


2The energy Vi? is proportional to C4 while accompanied by u?, so in the large p region this term 
may have some influence. 

3u, is the zero temperature correspondence of Ns, and it is not necessary to be a metastable vacuum. 
We also refer to it as Qs in the following. 


see it, consider the small deviation case and write us = (1 + ô)vs (|6| < 1), then the gap is 
approximated as 


A Viree = Vs View Va View 


2,2 
we h — CAA? v (u? — v?) + k =P 
1 
+zKAk [2u? (us — vs) + vs (v? — us) (3.12) 
i 
482K? (1+ 1/Rp) vt — 25 CAV? + T (3.13) 


Obviously, AViree goes to the doublet limit vm /4 as 6 + 0. In other words, a substantial 
deviation is necessary to decrease the energy gap away from the doublet limit. In fact, 
AViree can be even negative (we will see this soon later). A negative AViree is somewhat 
welcome since loop correction will be found to favor uplifting Vs relative to Vew. 

One could have a closer inspection into the deviation. Consider the minima structure 
of the singlet subspace at tree level, whose details are listed in Appendix. A, its absolute 
minimum locates at the origin or 


(1+ V1- 82x), (3.14) 


In Type-I the latter is just the case, which requires x, < 1/9. Using Eq. (3.8) one can 
rewrite x, as 


1 1 
Be = z — 3 (1+ Re) — Car’o?/ AQ, (3.15) 


In the limit CA — 0, one gets the following simple relation between us and vs: 


—v,(1+2/R,) +O0(C if Rk Z —1; 
won fl VEOGA 4 fa (3.16) 


Us +0O(Ca) if Rk S —1. 


which shows that u, usually deviates from v, significantly in the first case while in the 
second case they should be close to each other, given suppressed corrections from nonzero 
Ca. Note that at the leading order O(C4) is given by —C4A?v?/|1 + Ry|KAx, which 
indicates that the approximation breaks down for R,, near —1. In this case, a positive CA 
in xz, can also generate a deviation. 

Arguably, a substantial VEV deviation, i.e., for the first case in Eq. (3.16), tends to 
drive AViree < 0. Notice that in the decoupling limit C4 — 0, S = v, is always either 
a minimum or maximum (A, < 0 and R, Z —1) in the singlet subspace since the first 
derivative of V (S) from Eq. (3.10) over S$ is zero +. In the latter case, it is not surprising 
that Viiy(vs) > Vs(us); in the former case, a large negative Vs(us) is also possible for 
us S —vs (Rx > 0). Remind that the singlet-doublet mixing term is suppressed by small 
Ca, thus the above difference tends to dominate in AViree, rendering it negative. This is 
particularly true in the A, < 0 region when Vg(vs) is a maximum, where R,, > —1 requires 


“Turning on a small Ca will make Qpw stable and shift S = v, away from being a maximum. 


—10— 


—A,,/K > 4vs Z O(TeV) or even order of magnitude larger for a larger u. That large 
—A,,/K, in terms of the naive argument below Eq. (3.11), renders Qs well below OQpw. 
Therefore substantial loop corrections are indispensable to flip the order. 

Before heading towards the loop corrections, let’s make some observations of the tree- 
level results on the R, — AViree plane (see the upper panel of Fig. 2). They are in accord 
with the analysis above: (I) In Rẹ < —1 region, Us ~ vs, S0 AViree clearly takes the doublet 
limit vm? /4; (II) In R, ~ —1 region, the magnitude of AViree can blow up, in particular 
within the window —1 < R, < 0 and for a relatively large u; (III) In R, > 0 region, as 
argued before, AViree can also be negative and of order of a few 10° GeV“, significantly 
smaller than the case (II). More complementary analysis is left to the part of numerical 
study. 
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Figure 2. AViree and AVnum (with loop corrections) versus Rẹ. The left(right) two figures are the 
plots for the H,(H2)—scenario. Clearly, the large negative energy gap at the tree-level AViree is 
driven back to small positive value AVjum through loop corrections. 
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3.2.2 Loop-level analysis 


Previously it is shown at tree level Qg usually lies above Qgw. Here we will demonstrate 
that the tree-level order is going to be flipped by loop effects, which tend to lower Qgw but 
lift up Qg. In the DR scheme, the former is mainly ascribed to the remnant of Coleman- 
Weinberg potential after correcting the Higgs soft mass terms, while the latter is mainly 
due to the shift of us. In the following we describe their details respectively. 

On one hand, loop effects can lower Qgw. To offset the shift of VEV’s in Qgw due to 
Eq. (2.11), we may need to add the corresponding Higgs quadratic term, —5Am?¢?, with 
Am? determined to be 


p= 25 AO) neo], Liles) -2). (3.17) 


We have introduced A; = (—)?%n;/64n? and L; = ln mon for short. As a consequence, 
the remnant of Coleman-Weinberg potential, Vew — nee v7, results in a shift to the 
tree-level Vew vacuum energy: 


vw = Yo Amt (o) (Lio) i ley (Eale) —D]- (818) 


vı 
m} (vi) 
The above expression can be simplified greatly for two limits of m? (vi): 


Strong v;—dependence In this limit the mass of particle i dominantly originates in cou- 
pling to Higgs fields, such as the SM-particles and Higgsinoes. Then we have 


Vein (v1) =) es mi (v (z ilu) — 5) (3.19) 


Bare in mind that we have fixed Q = 2 TeV, thus the relatively light fermions, e.g., 
top quark and Higgsinos, contributes a positive Ve and make for flipping. The 
resulted decrease in Vgw can be up to order 10° GeV* for a heavy u ~ 500 GeV. 
By contrast, the light weak gauge bosons hamper for flipping but numerically it is 
unimportant for their lightness. 


Weak v;—dependence Some particles like stop have large (soft) mass terms, so they 
typically have quite weak dependence on v; °. In this limit one may write m? (vı) = 
m? (1+ f(x)) with z = v/m; and f(x) < 1. With that, we get an approximation 

; (Li — 3/2) — 

[(1 — L) (2f(x) — x f'(@)) — Lif(a) (f(x) - xf'(@))], (3-20) 


Aim; 


~a oh 


2 
Here L; = ln a The term in the first line is a constant thus contributing null to the 
energy gap. While for the second line, heavier CP even/odd Higgs bosons with their 


°In our paper we decouple the stop in the thermal plasma but include their radiative corrections to the 
Higgs potential. 
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mass dependences on the VEVs (f(x) ~ x) will benefit the reduction of Vé%,, at the 
VEVs. If f(x) = A?x? (stop without trilinear soft mixing), the leading z—dependence 


in the second row will vanish, with energy shift proportional only to powers of v4, 
i.e., Aj LX} v7. 


In summary, viewing from our particle spectrum, loop effects tend to decrease energy of 
Qgw. In the following we discuss the u,—shift effect on Vs. 

On the other hand, loop effects can lift Vs up. Here the discussion is different from the 
previous case, because in Qg the singlet VEV changes after loop corrections and the effective 
radiative potential plays an important role. It induces a shift to m, inherited from the 
previous discussions in Qgw. On top of that, it affects other tree-level couplings, as can be 
seen by expanding Vew(S) into polynomials of S. We give corresponding typical examples: 
The heavy Higgs bosons (still lighter than u) and Higgsinos with mass As increase me and 
k by an amount, respectively, ~ TeV?/167? and 


k (i Hiet i (3.21) 


8T2kK2 


From Eq. (3.14) and Eq. (3.11) we know that both |us| and |Vs| are monotonically decreas- 
ing functions of x, (and « as well from Eq. (3.15) in Ca — 0 limit.). Thus, in the region 
with relatively small x? (< 0.01), |us| may be decreased and negative Vg is increasing so 
Qs is lifted up. 

The loop-level numerical results on the Ry — AViree plane are shown in the lower panel 
of Fig.2. The second role of lifting Vs which effectively increases « is crucial to flip vacua 
order with an especially large tree-level gap, which is characterized by —1 < R, < 0. As 
mentioned before, it is usually accompanied with a relatively smaller « and large u, which 
yields a sizable increase of « from Eq. (3.21). Recall that R, « «, increasing x may drag 
Rp, out the window —1 < R, < 0 and make Rẹ < —1. Therefore from Eq. (3.12), It is 
reasonable for us to draw a conclusion, i.e., the loop-level gap goes to the doublet limit 
miv? /4 > 0 for those points. In this way, the tree-level order is flipped. 


3.2.3 Numerical results 

The SFOEWPT is the result of complicated interplay among quite a few parameters, 
including A, k, A, etc. We have turned to the numerical methods for a global understand- 
ing, and for a cross-check with previous qualitative analysis. Using the NUSSMTOOLS 
package [55-57], we scan the parameter space of the model with constraints from various 
relevant experiments, including the constraints on Higgs signatures. The parameter setting 
is listed as follows: 


k : (0.01, 0.5), A: (0.3, 0.8), tang: (1.5, 10), 
Ay : (200, 2000) GeV, A, : (1000, 1000) GeV, ju: (100, 600)GeV. (3.22) 


To minimize the uncertainty from the soft spectrum and explore the EWPT properties from 
the genuine Higgs-singlet sector, we assume that they are irrelevantly heavy. In particular, 
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the parameters in the stop sector are taken to be m2 = m2 = 2 TeV and A; = 0°. The 


numerical points used in the previous sections are actually abt this parameter space. We 
calculate the phase transition strength following the textbook way: Search for the minima 
of the complete one loop Higgs potential (2.10) at each temperature, and then find out 
Te and Ye by the degenerate vacua condition. In the following we show the parameter 


distributions favored by SFOEWPT, and try to give interpretations of them. 


H,-scenario Ho-scenario 


Yc 
o- NOK HOD 


-10-5 0 5 1015 20 25 30 35 -10-5 0 5 1015 20 25 30 35 
R R 


K K 


Figure 3. Rp versus ye in Type-I transition, with color code denoting u. Left panel: H;—scenario; 
Right panel: H2—sceanrio. 


Firstly, we present plots on the Rẹ — Yc plane in Fig. 3. From them one can see that Rx, 
is an helpful variable to judge ye. In other words, the previous tree-level analysis, despite 
of very complicated loop corrections, still provides valid insights. A possible large ye is 
accommodated in two distinct regions: 


Rs ~ O(1) — O(10) ye > 1 in Type-I favors the Rẹ > 0 (thus A,/« > 0) region, which 
could have a sizable (us — vs) deviation to avoid the doublet limit. One may wonder 
why large R,,’s fail. Consider R, >> 1 and still take C4 — 0, from Eq. (3.16) one 
gets us © —v,. Then it is straightforward to derive 


AVuree ~ MZV? /4 — 16K702/3R, + O(Ca). (3.23) 


So AViree, with its singlet contribution suppressed by the large Rx, goes to the doublet 
limit again. Numerical results show that in H2—scenario, ye > 1 occurs for Rẹ < 10, 
while in H;—scenario ye > 1 can still occurs for Ry, as large as 30. 


Note that R, can not be too small neither owing to phenomenological reasons. A 
large positive A, threatens the positivity of the light CP odd Higgs boson mass (we 
will turn back to it later). Moreover, a very small « leads to a light (singlino-like) 
neutralino, which opens too large 126 GeV Higgs exotic decay branching ratio in 


ĉIn fact, decoupling the stop is not merely requiring heaviness, and we additionally require a smaller 
A. This can be seen from Eq. (3.20), the presence of sizable A; would make heavy stops leave appreciable 
effects on energy gap at loop level. 
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the non-PQ limit. Besides, singlino-like LSP may be over-abundant because of too 
small annihilation cross section. In conclusion, A,/« can not be too large and ac- 
cordingly R, gets a lower bound, about 4.0 and 2.0 in H;—scenario and H2—scenario 
respectively as shown in Fig. 3. 


R, ~ —1 For R, < 0, it is not surprising that points with ye > 1 crowd around Rẹ ~ —1. 
On one hand, R, « —1 fails to achieve SFOEWPT owing to the doublet limit, as 
argued at the end of Section 3.2.1. On the other hand, R,, cuts off before approaching 
to zero because — A, is not allowed to be very large here owing to the singlet-like CP 
even Higgs boson. Therefore, R, is preferred to be around -1. 


Hı—scenario and Hj—scenario demonstrate a remarkable distinction in this region, 
namely ye 2 1 is accommodated in the latter but not in the former. , If the point is 
in the H,—scenario with ye Z 1, a large M33 term is required and from Eq. (2.8), we 
can see that this further requires a large My, (very large u and Ay, see Fig. 3 and 
Fig. 6) since the term 4«?v2(1 + 1/R,.) is close to its minimal. 


Notice that those points in H,—scenario have obviously large ws, which yield con- 
siderably loop corrections according to the discussion at the end of Section 3.2.2. 
Those loop corrections would bring ô back to zero (see Fig. 4) and therefore forbid 
the SFOEWPT. As a comparison, in H2—scenario ô could have a very large deviation 
from zero and even approach 6 > 1. It is precisely those points which have a small 
AV and trigger a SFOEWPT. 


H,-scenario H5-scenario 


0.1 


AV num(10°GeV) 
AV pum(10°GeV) 


0.01 
—2.5-2-1.5-1-0.5 0 0.5 


Ù ô 


Figure 4. Loop-level gap AVaum versus 6 = (us — Us)/vs for the Type-I transition for the Hı— 
and Hə—scenario respectively. 


Next, we outline the distributions of relevant parameters favored by ye > 1 in Type-I: 


e Light u and small tan 8 are preferred (see Fig. 5), in particular in H,—scenario where 
u S 250 GeV and tan 8 < 3.5. In H2—scenario they can extend to a bit larger regions, 
i.e., WS 250 GeV and tan 8 <5. In this sense, Type-I agrees with the most natural 
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NMSSM scenario [22], which have a slice of parameter space with A ~ 1, tanp ~ 1 
and u ~ mz. 


e A) ~ (200,600) GeV (see Fig. 6). Due to the suppressed mixing effect |C4| « 1, A) 
is strongly correlated with p, i.e., Ay © 2u/ sin 22. 


e As for k and A,, most of the preferences can be traced back to the discussion on 
R, = 4Kv;/A,.. We emphasize again that both H;— and H2—scenario accommodate 
SFOEWPT for A, > 0 , while for A, < 0, Ye can hardly achieve O(1) in H,—scenario 
(see Fig. 8). 


3.3 Type III: Results and analysis 


In this subsection, we turn our attention to the Type-III transition. This type of EWPT 
is studied in detail by the early works [45, 48-50] due to its compatibility with the near 


H,-scenario i H,—-scenario i 
600 600 
550 550 
500 500 
450 450 
400 400 
2 350 350 
300 300 
250 250 
200 200 
eo 150 > 150 
=— 100 100 
12345678910 1234567 8 910 
tanB tanB 
Figure 5. As in Fig. 3, plots on the tan 8 — ye plane. 
H,-scenario Hə-scenario i 
6 600 
550 
5 500 
450 
4 
400 
= 3 350 
300 
2 250 
1 200 
150 
0 0 100 
0 400 800 120016002000 0 400 800 120016002000 
Ay Ay 


Figure 6. As in Fig. 3, plots on the A) — Ye plane. 
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Figure 7. As in Fig. 3, plots on the « — ye plane. 
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Figure 8. As in Fig. 3, plots on the A, — Ye plane. 


PQ symmetry limit and its one-step nature. Here we revisit this type of transition in the 
spirit of energy gap. Most of the analysis is similar to that of Type-I. 

It is worth pointing out that Type-III arises not only in the case that the origin is 
indeed the absolute minimum in the singlet subspace but also in the case that the origin 
is metastable. The latter has a Type-I-like vacua structure but belongs to Type-III due 
to thermal evolution: At T = 0, the absolute minimum in the subspace locates at us 4 0, 
but as temperature increases it will exceed the origin and become energetically disfavored. 
So the degenerating eventually happens between the origin and Qpw, who determine the 
gap. Such type-crossing phenomenons are not difficulty to be understood. The critical case 
Tk = M/A? = 1/9 < 1 means that m2 is positive and small (typically ms < 100 GeV), 
and thus sensitive to temperature. Although it is unable to distinguish quantitatively this 
case from Type-I at T = 0, our numerical results tell that the models of this case have a 
smaller gap between the origin and Qg than those in Type-I. 
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We proceed to discuss the energy gap at tree and loop-level. This time, the tree-level 
gap is quite simple, given by —Vew. 


A Vinge => View Viti View 
vii +C eee = (3.24) 
4 ‘Apu + 620% 3R, : 


Obviously, the singlet part dominates the energy gap for a large u, and thus a small gap 
appeals to a large negative A, such that —4/3 S R, < 0. For a moderate u, the mixing 
part becomes important and a negative C'4 can help to decrease the gap, which requires a 
large Ay. The features outlined above are well reflected in the Fig. 9. 

The above analysis can be well adjusted in the near PQ-limit, where the mixing part 
from the second term in Eq. (3.24) could play a dominate role, so a large A, and a 
moderate u play key roles in decreasing the energy gap. Notice that in Ref. [50] where a 
relatively larger tan 8 > 10 region is considered, the corresponding A) ~ (2,5) TeV is even 
larger in order to achieve a small AV for SFOEWPT. We observe that Type-III can only 
be accommodated in a very restricted region. This is because that as soon as R, become 
a bit larger, Qs will become so deep that the transition changes into Type-I ”. 

The loop-level analysis in Section 3.2.2 is also applicable to Type-III, i.e. there are 
mainly two kinds of loop-level effects. On the one hand, loop corrections lower Qgw and 
make it to be the global minimum in the same way as in Type-I. On the other hand, the 
loop corrections in the singlet subspace lift up Qs by increasing m and «. However, the 
latter has no influence on the energy gap concerned here, and just makes Qs shallow and 
ready to exceed the origin in thermal evolution. 


H,-scenario H»-scenario 


2.5 -2 -15 1 -05 0 5 4 -3 -2 -1 0 


Figure 9. Rẹ versus yc in Type-III transition, with color code denoting u. Left panel: H;—scenario; 
Right panel: H2—sceanrio. 


Finally, we summarize the parameters preference of SFOEWPT in Type-II: 


“the points in R, < —1 region have a negative C4, which can relax the conflict between large Rẹ and 
Type-III, as indicated by Eq. (3.8) or Eq. (3.15). 
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e Unlike Type-I, SFOEWPT in Type-Ill prefers a larger u = 250 GeV, especially in 
the region —4/3 < Ry, < 0 (see Fig. 9). As argued above, AV here is effectively 
decreasing by the singlet part. Similar to Type-I, tan 8 < 3.5 is favored (see Fig. 10). 


e The favored values of Ay here is larger than in Type-I since A) ~ 2y/ sin 28 or even 
larger for a negative C4(See Fig. 11). 


e Another obvious distinction between Type-I and -III can be observed in Fig. 13: 
A,, > 0 barely accommodates Type-III, not mentioning to realize ye > 1. Eq. (3.8) 
may provide a simple interpretation. It indicates that given a positive A, we have 


2 aie 


ine => —KA,Us — Kv < 0 and consequently x, < 0, which yields a large deep 


minimum in the singlet direction, therefore strongly favors Type-I. 


e The vast majority of surviving points are in Hy—scenario. As a matter of fact, 
seemingly H,—scenario has a strong tension with Type-III. This may blame to its 
parameter configurations that is difficult to achieve a large M33 unless u and A) are 
large, which forbids the SFOEWPT just like the case of Type-I. 


Ye 


H,-scenario i Ho-scenario i 

600 600 

550 550 

500 500 

400 400 

350 350 

300 

si 250 

250 200 

200 150 

= 150 100 

12345 67 8 9.10 12345678910 
tanB tanB 


Figure 10. As in Fig. 9, plots on the tan 8 — ye plane. 


4 Relevant Phenomenologies 


4.1 Higgs Spectra 


The most remarkable connection may lie in the existence of a relatively light Higgs boson 


in the spectrum, though their masses are heavier than mp,,,/2 otherwise the Higgs exotic 


decay channel would dominate and get severely constrained from the current Higgs signal 


rates. 


For the CP odd Higgses, it is well known that the absence of tachyon states in the 


Higgs CP odd sector strongly favors negative kA, and now it is more constrained from the 


observed Higgs signal rate as mentioned above. Such a sharp contradiction substantially 
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Figure 11. As in Fig. 9, plots on the A, — Ye plane. 
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Figure 12. As in Fig. 9, plots on the «K — ye plane. 


compresses the allowed parameter space and leads to a remarkable prediction on the Higgs 
spectrum. The CP odd Higgs sector contains two physical states a12. In the case of 
M? > M2 we have 

ma, ~ A (Ay + 4kus) Vuva/Us — 3K ARVs (4.1) 


where the second term can be expressed as —12«7v?/R,. Thus, a positive small Rẹ drives 
the lighter CP odd Higgs boson mass downwards. 

For the CP even Higgses, R, can also affect the singlet-dominated CP even Higgs mass, 
as we mentioned in Eq. (2.8): 


1 
(M2)as = 4n202 (1 + x) iss (4.2) 
Rpg 
where --- stands for other terms either are small or not relevant here. Therefore the 


preference of R, from SFOEWPT may result in some specific distribution in the Higgs 
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Figure 13. As in Fig. 9, plots on the A, — Ye plane. 


spectrum. Especially we can see that —1 < Rẹ < 0 is the region that the above term 
get its minimal and will drive down the hı mass significantly. Therefore, it is clear that 
in the NMSSM, the SFOEWPT will impose a specific Higgs spectra through the critical 
parameter Rx: 


e R, > 0 prefers a light CP odd Higgs mass (small m1) with no strong preference on 
Mhl. 


e 0 > R, > —1 prefers a light CP even Higgs mass (small mp; in the Hə—scenario) 
with no preference on mq. 


The above speculation is confirmed by the histograms Fig. 14 and 15 in which we have 
imposed the CMS Higgs signal data on surviving samples in various channels including 
hsm > Z2Z2,WW, yy, bb,T7 at 2o level [58] (ATALS constraints have a relative smaller 
number of points with the same distributions). In Fig. 14, we compare the normal his- 
tograms for Mp1, Ma, and those with SFOEWPT. We can see that the latter has more 
points concentrated on the light mp1, mq, region. For those histograms which we have 
both R, > 0 and Rẹ < 0, we further distinguish them in the histograms Fig. 15. We can 
see that there is a very distinction that R, > 0 prefers small maı while Rẹ < 0 prefers 
small mp1, which is coincident with our analysis above. 


4.2 Dark Matter Consideration 


If the lightest NMSSM neutralino is assumed to be the WIMP Dark Matter (DM) candi- 
date, an important consideration is the DM relic density in the present epoch. Combining 
the PLANCK [59] and WMAP 9-year data [60] and also including a 10% theoretical un- 
certainty, the 20 range of the WIMP DM relic density can be considered in the following 


range: 


0.091 < Qh? < 0.138 (4.3) 
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Figure 14. Left panels: Histograms of mai for physical points satisfying the CMS Higgs signal 
at 2c level (blue-blank) and SFOEWPT (red-shaded) for Type-I samples in H2—scenario (upper) 
and H,—scenario (bottom) . Right panel: As in the left panels, histogram of mg, in H2—scenario 


is plotted. 


Since the bino and wino mass parameter M1, Mə have been fixed at 2 TeV in our analysis, 
the lightest neutralino basically consists of higgsino and/or singlino. From the neutralino 
mass matrix one can tell that, « can significantly affect the singlino component in the DM: 


g _ h 
Mo n T 0 
Mo = 0 =u —Avu (4.4) 
0 Avg 
2KUs 


Since all squarks and sleptons have decoupled in our analysis, the DM can annihilate in 
the early universe only through the light Higgs bosons in the s-channel process or as the 
final states when the process is kinematically opened. For a highly higgsino-like DM, the 
coupling of the Higgs with the DM can be sizable and it is very easy to obtain a small 
relic density. We have checked that for each SFOEWPT scenario in our discussion, it is 


— 
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Figure 15. Left panel: Histogram of Mmaı for physical SFOEWPT points with arbitrary R,.’s 
(blue-blank) and points with R, < 0 (red-shaded) for Type-I transition in H2—scenario. Right 
panel: As in the left panel, histograms of m, 1 is plotted. 


always possible to pick out several samples which can produce a relic density that does not 
overclose the universe, or even lies in the band shown in Eq. (4.3). 


5 Conclusion and discussion 


After the discovery of 126 GeV Higgs boson, the NMSSM is an attractive supersymmetric 
theory in virtue of its specific tree-level effect to enhance Higgs boson mass and allowing 
a more natural u parameter. On top of that, the tree-level effect can readily enhance the 
strength of EWPT ye 21.0, which is required for a successful EWGB mechanism to generate 
the baryon asymmetry. In this article we have concentrated on studying SFOEWPT in 
the NMSSM, paying special attention on its relation with Higgs phenomenology. We have 
calculated EWPT strength ye with the one-loop finite temperature effective potential and 
find that a larger ye requires a smaller gap AV. Then, in terms of the vacua structure 
and its evolution with temperature, we divide EWPT into three categories, Type-I, II and 
II along with two Higgs spectra patterns: H;—scenario and H2—scenario. We use our 
semi-analytical analysis as the intuitive understandings and then use our numerical results 
to confirm those understandings. 

We have observed a dimensionless critical parameter R, = 4Kvs/A, which has demon- 
strated a clear correlation between the different types of SFOEWPT in NMSSM and the 
Higgs spectra as follows: 


e In Hı—scenario, the Type-I phase transition prefers R, > 0(C (5,30)) and much few 
points exists in the Type-III phase transition with small negative R,,, 


e In Hə— scenario, the Type-I phase transition has two distinct regions, in which either 
0>R, > —1 or R, > 0(C (2,7)). For the Type-III phase transition, much of them 
lies in the region 0 > Rẹ > —4/3. 
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A SFOEWPT in general prefers a relatively light CP odd or even (H2—scenario) Higgs 
mass. In particular, R, > 0 prefers small may while Rẹ < 0 prefers small mpy. 

The current classification on the EWPT patterns and the Higgs spectra has a great 
importance to guild us to the next step of understanding the nature of EWPT. In particular, 
the Higgs spectra in SFOEWPT in the NMSSM, or even broadly, in some SM / 2HDM + 
singlet models, prefers either a light CP odd or even Higgs with their mass slightly more 
than mp,,,/2. We hope that our work can help bring more attention on the observation 
of those (60, 100) GeV light Higgs and works along this direction will be presented in the 
future 8. 
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The range of parameters minimum 
I 

t >w > t 0 
by > 0 1 b3 <0 Sa 

< 
pal 7Sa Ts>0] 5 
b3 <0 Sa 

b2 < 0 

° bs > 0 S 


Table 1. Minimum structures in singlet direction 


A Minimum in the singlet subspace 


For convenience, we present the minimum structure in the singlet subspace here. The 
potential has a form 
b3 


— b2 Go 3, b4 
=a rae Tr] 


V(S) st. (A.1) 


The minimum condition 0OV/0S = 0 has following solutions: 


S=0 (A.2) 
—be + 2 
a va Ababa (A.3) 
4 


The latter is physical as long as A = b — 4bgb4 > 0, i.e, £ = b2b4/(2b2) < 1/8. 
The minimum structures are listed in Table 1. Sometimes, it is useful to express the 
potential at nontrivial minimum as 


Obh S 


2 — 
-E+ -ga O) 


Vinin ( St ) 


(A.A) 


where f(x) = 4+ 4/1 — 82 — 482 — 322\/1 — 82 + 9627, which decreases monotonically 
from one to zero as x increares from zero to 2/9. Obviously, Vmin = 0 for x > 2/9. 
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